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DERIVATION AND TESTS OF THE GODDARD COMBINED 
GEOPOTENTIAL FIELD (GSFC 1.70-C) 


INTRODUCTION 

The potential of the earth at a point with spherical coordinates (r,0,k), in an earth centered 
rotating coordinate system may be given by 



co ( a ^ "| 

£ £(7) p t . <sin Q (q e. » cos m x + s -t, » sin m X) r 


(i) 


f.= 2 m " 0 


This form of the potential is the recommended form for axially asymmetric (zonal, tesseral, and 
sectorial harmonics) cases, Reference 1. In equation (1), fx is the product of the gravitational 
constant, G, times the mass of the earth, M, and a e is the mean equatorial radius of the earth. 
The Legendre associated function, m (sin<£ ), is defined by 


? 4, m = 


2 l l \ 


a-* 2 ) 


2 \tn /2 


(X 2 - 1) 


l 


d x 


'&+ 1 


( 2 ) 


Instead of using the spherical harmonic coefficients C^ and S|> ; the fully normalized coef- 
ficients C^ m and have been adopted for this work. They are related to C^ and coef- 

ficients by 


( 3 ) 


= Np C. 

As, m At, m As, 


where 




So - No So 

As, m As, m At, 


K (2 l + 1) (l - m) ! 


(f- + m) 1 


- 1/2 


and where K = 1 for m = 0, and K = 2 for m d; 0. 


Combined geopotential fields based upon four individual solutions for the geopotential have 
been published, Reference 2. The first such field, C, appearing in Reference 2 was obtained by 
taking an arithmetic mean of these four solutions which were obtained from satellite tracking data. 
The second field, CA, was obtained by taking field C and combining it with terrestrial gravity data. 
In this paper Professor Kaula stated that an arithmetic mean of different solutions from satellite 
data is superior to any single solution and a combination of satellite and terrestrial solutions, 
such as solution CA, should be superior to either solution alone. He also suggested that such a 
combination should be made when more solutions are available. 


In this paper we will form such a combination solution based upon eleven complete solutions 
and four partial solutions for the geopotential. This solution is not completely analogous to either 
of the solutions C or CA of Kaula. It is an arithmetic mean of individual solutions as is solution 
C; however, the individual solutions are made up of determinations based upon satellite data and 
ones based upon both satellite and terrestrial data. Thus, this combined solution is an arithmetic 
mean of solutions of the type C and CA. 

INPUT GEOPOTENTIAL FIELDS 

The fields used in obtaining the combined geopotential field are of various types. They include 
fields obtained solely from satellite data. These satellite determined' fields include ones obtained 
from optical observations, Doppler observations, and both optical and Doppler observations. There 
are also fields represented which were determined from both satellite and terrestrial gravity data. 
In addition to these two types of fields, various sets of coefficients determined from satellites in 
resonant orbits were used in obtaining the combined field. The particulars concerning these input 
fields are contained in Table 1. 


Table 1 

Gravity Fields used in the Combination Solution 



Derivation of Fields 


Gravity Field 

Reference 

No. of Satellites 

Data Types 

Comments 

NWL 5E-6 

4 

3 

Navy Doppler 

Complete to (7,6) 

APL 3.5 

5 

5 

Navy Doppler 

Complete to (8,8) 

SAO M-l 

3 

16 

Optical 

Complete to (8,8) 

Kob'dein 

6 

16 

Optical & Terrestrial 

Complete to (15,15) 

Kai a UCLA 656 

7 

9 

Navy Doppler, Optical 
Si Terrestrial 

Complete to (8,8) 

Kaula CA 

2 


Navy Doppler, Optical 
St Terrestrial 

Complete to (7,2) 

Kaula K-8 

2 



Complete to (7,5) 

SAO B6.1 

8 

24 

Optical, Range 
& Range Rate 

Complete to (16,16) 

SAO B13.1 

9 

■ 

24 

Optical, Range, Range 
Rate Si Terrestrial 

Complete to (16,16) 

Rapp 1967 

10 

16 

Optical Si Terrestrial 

Complete to (14,14) 

Rapp 1968 

11 

16 

Optical Si Terrestrial 

Complete to (14,14) 

Wagner 

12 

3 

Mean Kepler Elements 

Resonant Coefficients (2,2), (3,3) 

Gaposchkin St Veis 

13 

2 

Optical 

Resonant Coefficients (13,12), (14,12), (15,12) 

APL 

14 

3 

Navy Doppler 

Resonant Coefficients (13,13), (15,13), (17,13) 

Goddard 

15 

■ . . ' . 

. 

1 

Optical Si Range 

Resonant Coefficient (14,13) to be used with APL 
coefficients in April 1968 for GEOS II 


NASA.GSFCT&OS 

MISSION & TRAJECTORY ANALYSIS DIVISION 
BRANCH 552 DATE 04/03/70 
B Y Murphy & Marsh PLOT NO. 2099 ’ 

The result of taking the arithmetic mean of these fields appears in Table 2. Note however, 
that the combined field, GSFC 1.70-C, so obtained was truncated after the terms of fifteenth degree 
and order. If terms of higher degree were averaged, there would be only a few contributors. A 
variety of tests, conducted using this combined field are discussed in the next section. Tests were 
made with the remaining geodetic constants (zonals, earth radius, etc.) equal to those adopted for 
the SAO 1966 Standard Earth, Reference 3. For convenience they are also listed in Table 2. 
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Table 2 
Earth Model 

TCSSERAL HARMONICS* 


N 

2 

M 

2 

C< N, M) 
2.402 

S (N.W ) 
-1.370 

3 

1 

1 .891 

0 .219 

3 

2 

0.334 

-0.653 

3 

3 

0.5 18 

1 .355 

4 

1 

-0.550 

-0.452 

A 

2 

0.316 

0.604 

4 

3 

0.915 

-0.121 

4 

4 

-0.125 

0.212 

5 

1 

-0.037 

-0.061 

5 

2 

0.553 

— 0 o 24 8 

5 

3 

-0.357 

-0.012 

5 

4 

-0.142 

0.117 

5 

5 

0.040 

-0.501 

6 

l 

- 0.086 

0.055 

6 

2 

0.027 

-0.329 

6 

3 

0.051 

0.079 

6 

4 

-0.090 

-0.459 

6 

5 

-0.220 

-0.50 l 

6 

6 

-0.072 

-0.261 

— 7 

1 — 

U. 1 59 

0.037 

7 

2 

0.324 

0.089 

7 

3 

0.206 

-0.094 

7 

4 

-0.224 

-0.043 

7 

5 

0.055 

-0.040 

7 

6 

-0.266 

0.115 

7 

7 

o.o >-0 

0 70 43 

3 

1 

-0.047 

0.023 

3 

2 

0.059 

0.020 

8 

3 

-0.012 

0.081 

8 

4 

-0.094 

0.025 

8 

5 

-0.050 

0.02 1 

5 

6 

— U • u 3o 

0.259 

8 

7 

0.040 

0.025 

3 

8 

-0.142 

0.020 

9 

1 

0.111 

-0.014 

9 

2 

0.0 03 

-0.008 

9 

3 

-0.034 

-0.030 

~9 

4 

0. U 6' 0 

0.053 

9 

5 

-0.044 

-0.054 

9 

6 

0.062 

0.035 

9 

7 

-0.055 

-0.087 

9 

8 

0.218 

0.037 

9 

9 

0.012 

-0.025 

TO 

l 

o * (j rtj 

— U n 0 t> 4 

10 

2 

-0.067 

-0.034 

10 

3 

-0.017 

-0.063 

IQ 

4 

-0.050 

-0.125 

10 

5 

-0.045 

0.006 

10 

5 

-0.076 

-0.060 

10 

7 

0. 07 T 

~ : — 0.013 

10 

8 

0.079 

-0.097 

10 

9 

0.040 

-0.044 

10 

1 0 

0.051 

-0.064 

1 1 

1 

-0.034 

-0.027 

11 

2 

0.057 

-0.033 

1 1 

3 

-0.023 

-0. 103 

1 1 

4 

— Q • 0 0 0 

-0.019 

1 1 

5 

0.041 

-0.006 


*UN1T EQUALS TEN TO 
MINUS SIXTH POWER 


N 

M 

C f N ♦ M ) 


l 1 

6 

-0.044 

-0.073 

1 1 

7 

0.012 

-0.101 

1 1 

8 

0.160 

0.014 

1 1 

“ 9 

- 0 . 0 l b 

"-ir.im' 

1 1 

10 

-0.046 

-0.040 

1 1 

11 

0.0 83 

0.012 

12 

1 

-0.079 

-0.047 

12 

2 

0.004 

0.046 

12 

3 

0.073 

-0.01 1 

12 

4 

- 0 . U 39 

— u • uy b 

12 

5 

—0.002 

0.065 

12 

6 

-0.065 

0.024 

12 

7 

0.004 

0.026 

12 

8 

0.009 

0.073 

12 

9 

-0.073 

0*0*3 

12 

1 U 

— OTtJTTZ” 

y.oO'j 

12 

1 1 

-0.029 

-0.011 

12 

12 

-0.025 

-0.012 

13 

. 1 

-0.017 

0.01 1 

13 

2 

-0.031 

-0.031 

13 

3 

-0.067 

0.024 

1 •> 

4 

—0.025 

0.039 

13 

5 

0.083 

-0.077 

X3 

6 

-0.038 

0.654 

13 

7 

0.016 

-0.007 

13 

8 

0.047 

-0 .025 

13 

9 

-0.010 

0.009 

L 3 

1 0 

0.055' 

— 0.07 1 

13 

1 1 

-0.033 

0.009 

13 

12 

-0.030 

0.072 

13 

13 

-0.093 

0.055 

14 

1 

-0.0 17 

0.044 

14 

2 

-0.028 

0.009 

1 4 

3 

0.109 

— 0 . 0 2 1 

14 

4 

0.032 

-0.043 

14 

S 

0.068 

“0.039 

14 

6 

-0.033 

-0.01 1 

14 

7 

0.051 

-0.015 

14 

8 

-0.032 

-0.056 

1"4 

" " 9' 

-0 .027' 

0.054 

14 

1 0 

0.009 

-0.070 

14 

1 1 

0.051 

0.02* 

14 

1 2 

0.017 

-0.027 

14 

1 3 

0.023 

0.051 

14 

1 4 

-0.040 

-0.009 

1 5 

l 

— 0.0 1 5 

-0.037 

15 

2 

0.004 

-0.109 

15 

3 

-0.033 

0.030 

l 5 

4 

0.023 

0.074 

15 

5 

0.022 

0.038 

15 

6 

0.050 

-0.053 

T5" 

. y 

0.075 

0.035 

15 

0 

-0.030 

0.031 

15 

9 

-0.040 

0.045 

15 

1 0 

0.038 

0.022 

15 

1 1 

0.0 13 

0.072 

15 

12 

-0.048 

0.032 

' l 5 ; 

1 3 

*‘“0.0 43 

'—"O's "Q'33 

15 

14 

0.0 08 

-0.013 

is 

15 

-0.005 

-0.013 


Z3NAL HARMON ICS*C< N« o) 


C( 2.0)=- 

43^.1733 

C( 3.0)= 

0.9623 

C( 4.0)= 

0.5497 

C ( 5.0)= 

0.0633 

C( 6,0)= 

-0 • 1 79? 

C ( 7.0)= 

0 .0360 

C( 8,0)= 

0.0655 

Z( 9,0)= 

0.0122 

C( 10 * o ) = 

0.01 13 

fill . 6) = 

076636 

C ( 1 2 * 0 ) = 

0.0714 

CC l 3 » 0 ) = 

0.0219 

C ( 1 4 . 0 ) = 

-0.0332 

C ( 1 5 . 0 ) = 

0.0 


SEMI MAJOR AXIS 
(METERS) 

6378U2.00 


FLATTENI NG 


1 . /29S.25 


GRAVITATIONAL CONSTANT 
( METERS* *3/ SECONDS ) 

3 • 986009000 14 
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TESTS AND EVALUATION OF THE COMBINED FIELD 


Several tests of this new model have been made. The first pair of these were concerned with 
the fit to six day arcs of precision reduced optical data for GEOS I and GEOS H. The fits were 
made without and then with special values for resonant coefficients for the two satellites. For GEOS I 
the resonant coefficients were those obtained by Gaposchkin and Veis, Reference 13; and for 
GEOS II the resonant coefficients were those obtained by APL, Reference 14, and Goddard, Refer- 
ence 15. Characteristics of these two satellite orbits are given in Table 3. In Table 4 there appears 
a list of resonant coefficients obtained previously for these two satellites. The combined field was 
then tested in orbit computations with and without resonant terms against some of the fields that 
were used in deriving it. After these tests were completed, the combined field was truncated in 
the following manner. All the coefficients of degree greater than eight with the exception of those 
of order twelve, thirteen, fourteen, and fifteen were deleted from the field. The same orbit deter- 
mination runs were then repeated for this adjusted field, GSFC 1.70-T. 


Table 3 

GEOS I and GEOS II Orbit Characteristics 



GEOS 1 

GEOS II 

semi major axis 

8077.9 km 

7701.1 km 

eccentricity 

.070 

.033 

inclination 

59?4 

105?8 

argument of perigee 

312 ?0 

194°8 

right ascension of 

266?1 

35377 

ascending node 
mean anomaly 

224?1 

121?0 

epoch 

July 11, 1966 0 h 0 ra 

April 28, 1968 17 1 ' 56 m 

resonant harmonics 

12th order 

13th order 

resonant period 

7.3 days 

6.5 days 

perigee height 

1130.6 km 

1071.8 km 

apogee height 

2268.9 km 

1574.1 km 
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Table 4 

Resonant Coefficients for GEOS Satellites 


GEOS I* 



GEOS II** ] 

Ci2, 12 

= 

-3.1 X 10" 8 

G 13, 13 


-6.53 x 10 _a 

S 12, 12 

= 

.08 x 10 -8 

S,3.13 


5.79 X 10“ 8 

C 13,12 

= 

-6.769 X 10" 8 

G 15, 13 

= 

-3.82 X 10" 8 

S 13, 12 

= 

6.245 x 10~ fi 

S 15,13 

.= 

-1.85 X10" 8 

C 14 , 12 

= 

.261 x icr 8 

G 17, 13 

= 

6.36 X 10~ 9 

S 14 , l 2 

= 

-2.457 X 10“ 8 

S 17, 13 

= 

1.11 XI O' 8 

C 15, 12 

= 

-7.473 X 10“ 8 

C 14,13 

- 

7.81 X 10" 9 

S 15,12 

= 

-1.026 x 10" 8 

s 14, 13 


8.91 X 10“ 8 
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SAO 

"APL & GSFC 




The fits to the data for the six day GEOS I and GEOS II arcs with the GSFC 1.70 C and GSFC 
1.70T geopotential fields appear in Table 5. Also presented in Table 5 are fits obtained from 
previous studies, Reference 16, for the same orbital arcs using some of the complete fields used to 
derive the combined field. For the orbit computations involving the SAO B13.1 field, the recom- 
mended set of zonals, Reference 17, adopted for the 1969 SAO Standard Earth was used. Two 
points can be made concerning this table. First, after the resonant coefficients for these two satel- 
lites are inserted into the field, the fit to the data with the combined field is better than that of any 
other field for GEOS I, and for GEOS II the fit using the combined field was second only to the SAO 
B13.1 field. Secondly, after the field is adjusted so that it has only about forty percent of the terms 
that the original field had, the fit to the data remains about the same. 


Table 5 

Rms’s (Sec’s of Arc) of Fitted Orbits for Six Day Arcs of Optical 
Data from GEOS I and GEOS II when Different Gravity Models 
are used with and without Resonant Coefficients* 


Geopotential 

Model 

GEOS I 

GEOS II 

Number of 
Coefficients 

without 

with 

without 

with 

SAO M 1 

19.04 

2.52 

17.36 

3.08 

108 

SAO B13.1 

2.56 

3.27 

3.39 

2.29 

294 

Kbhnlein 

14.65 

2.89 

9.41 

3.12 

236 

Rapp 1967 

7.81 

6.91 

11.30 

5.48 

206 

NWL 5E-6 

11.82 

3.33 

27.99 

8.08 

58 

APL 3.5 

13.51 

6.64 

59.41 

5.79 

74 

Kaula 1967 

5.80 

5.95 

16.67 

9.32 

88 

GSFC 1.70C 

6.48 

2.30 

6.67 

2„95 

236 

GSFC 1.70T 

6.67 

2.43 

6.87 

3.23 

88 


NOTE 1" equals approximately 7 m. 
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Certain coefficients in the combined field were adjusted using optical data in the six day arcs 
referred to in Table 5. Other coefficients were adjusted using accelerations from seven twenty- 
four hour satellites in the Syncom, Intelsat and ATS (Applications Technology Satellite) series of 
spacecraft. These accelerations had been prepared by Mr. C. Wagner for use in some, of his deep 
resonance studies (see Reference 12, for example) and made available for this work by him. Values 
of the original and adjusted geopotential coefficients using these data appears in Table 6. The 
second and third degree sectorial harmonics were improved with the sectorial harmonic of degree 
four and the tesseral harmonics of degree three order one and degree four order two which are 
also sensitive for the accelerations held fixed to the values for the coefficients appearing in the 
combined field. 


5 



Table 6 


Adjusted Coefficients* 


Term 

A Priori 

Improved 

Satellite (S) 

C2 , 2 

2.402 

2.434 

ATS 1,3,5, INTELSAT 1,2-F3, Syncom 2, 3 

S2.2 

-1.370 

-1.398 

* 

C3 , 3 

0.618 

0.726 


® 3 ,3 

1.365 

1.371 


^ 13, 12 

-0.0675 

-0.0703 

GEOS I 

^ 13 . 12 

0.0622 

0.0681 


e 

J 3,13 

-0.0653 

-0.0629 

geos n 

- 13, 13 

0.0579 

0.0586 

.... . 
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When the twelfth order coefficient was adjusted, a two second of arc fit was obtained for the 
GEOS I arc. The corresponding improvement to the fit for the GEOS II arc was only a few percent. 
The RMS error to the twenty-four hour accelerations for the low degree coefficients that are sensi- 
tive to these accelerations appears in Table 7 for the SAO 1966 Standard Earth, SAO B13.1, SAG 
1969 Standard Earth, GSFC 1.70-C and GSFC 1.70-C with the adjusted sectorial harmonics. 

It may be concluded that based upon the results so far, the GSFC 1.70-C field modified with 
the SAO, APL, and Goddard GEOS I and If. coefficients, Table 4, further modified with the adjusted 
coefficients of Table 7 might be the best field to be presented here for satellite orbit computations. 
If speed of computation becomes a factor, the analogue of this field based upon GSFC 1.70-T might 
be used with some small sacrifice in accuracy of computations instead of GSFC 1.70-C, 


Table 7 

Fits to 24 hr Satellite Accelerations 


Geopotential Field 


1966 SAO Standard Earth 

8.5 x lO' 7 

SAO B13.1 

3.9 x 10" 7 

1969 SAO Standard Earth 

2.5 X10" 7 

GSFC 1.70-C 

3.8 x 10' 7 

GSFC 1.70-C + Adjusted Sectorials 

1.9 x 10" 7 
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As a means of further testing the orbital solutions obtained with the GSFC 1.70-C and the 
modified GSFC 1.70-C gravity models, satellite position differences were computed using the 
orbits based upon the SAO M-l (1966) model (modified by GEOS-I and GEOS-II t esonant terms) 
and the more recent SAO B13.1 model as standards. 

The position differences were computed at five minute intervals and were resolved into radial, 
along track, and cross track components. The differences were computed using the unit vectors 
H, L and C j respectively which were calculated from the following relationships: 


R 

H = — 
- R 



C = HxL 

is the vector from the geocenter to the satellite, 
is the distance from the geocenter to the satellite, and 
is the velocity vector of the satellite. 

The RMS of satellite position differences for the respective GEOS-I and GEOS- II orbital arcs 
are presented in Table 8. The largest position differences were along track for both the GEOS-I 
and GEOS-II orbits which was not unexpected due to poorly modeled resonance. These along track 
position differences were reduced significantly when the GSFC 1.7G-*C model was modified by the 
Gaposchkin-Veis 12th order coefficients for GEOS-I and the APL and Goddard 13th order co- 
efficients for GEOS-II. As indicated in Table 8, the RMS of the along track position differences 
between the GEOS I and GEOS II SAO M-l standard orbits and the orbits computed with the modi- 
fied GSFC 1.70-C model were on the order of 25 meters. The RMS of the radial and crosstrack 
position differences were on the order of 10 meters or less for all cases. The comparisons be- 
tween the orbits computed using the GSFC 1.70-C modified models with the coefficients (13,12) 
adjusted for GEOS-I and (13,13) adjusted for GEOS-II showed an increase in the along\track com- 
ponent of 8 meters for GEOS-I and 1 meter for GEOS-II. This is consistent with the RMS of fits 
presented in Table 5 and the fact that adjustment of (13,13) for Geos II resulted in very little change 
to the orbital fit. 

Although the RMS of fits for the GEOS-I orbital arc were within .3 arc seconds for the SAO B13.1 
orbit and the GSFC 1.70-C orbit (2.56 arc seconds vs. 2.30 arc seconds as shown in Table 5) the 
trajectory comparison presented in Table 8 shows that the RMS of position differences were on the 
order of 67 meters. 


where 


R 

R 

v 


Ta.ble 8 

Satellite Position Comparisons 
GEOS I July 11“16 ; 1966 
GEOS II April 28-Ma.y 4, 1968 


Satellite 

SAO M 1 (Modified)* 

RMS Position Difference (meters) 

Maximum 

' 

VS 

Radial 

Crosstrack 

Alongtrack 

Total 

Difference 

GEOS I 

GSFC 1.70-C 

7.0 

7.8 

82.8 

83.4 

189.7 

GEOS I 

GSFC 1.70- C + 







Gaposchkin and Veis 12th order terms 

9.0 

6.4 

25.8 

28.1 

71.5 

GEOS II 

GSFC 1.70-C 

10.3 

9.3 

91.8 

92.9 

210.8 

GEOS n 

GSFC 1.70-C + 

7.8 

8.3 

26.9 

29.2 

75.1 


APL (1968) and GSFC (1969) 13th order terms 

. 


* 




SAO B13.1 







vs 



. 


' _ 

GEQS I 

GSFC 1.70-C + 

19.9 

7.3 

59.0 

62.7 

159.3 


Gaposchkin and Veis 12th order terms 






GEOS IT 

GSFC 1.70-C + 







■ APL (1968) and GSFC (1969) 13th order terms 

20.3 



14.3 

75.6 

'78.6 

188.2 


* 

Modified with Gaposchkin and Veis 12th order terms for GeOs I Comparisons. 
Modified with APL (1968) 13th order terms for Geos 1! Comparisons. 
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To this point the GSFC 1.70C field has been tested against other fields using several satellites 
in distant 24hr orbits and the GEOS satellites in their close, resonant and drag free orbits. The 
next test of this model will involve orbit computations using tracking data from the Orbiting Geo- 
physical Observatory (OGO-4) satellite. Fits to six two day arcs of Minitrack and Goddard Range 
and Range Rate data using NWL 5E-6, SAG Ml, SAO B6.1, and GSFC 1,.70C are presented in Table 9 

The. analysis of OGO-4 data is of particular interest for a number of reasons. In the first 
place OGO-4 was not used in deriving any of the models discussed in this paper. This satellite 
in its near polar close earth orbit (a=7023.4km, e=,03492, 1=86? 00) with only modest drag and 
slight resonance (150 m) is ideal for sampling the total gravity field. 

Orbital comparisons similar to the ones made for the GEOS orbits have been also made using 
the OGO-4 arcs discussed above. In this case We were interested in the overlap region of succes- 
sive pairs of orbital arcs. The results of these orbital comparisons appear in Table 10 for the 
NWL 5E-8, SAO Ml, SAO B6.1 and GSFC 1. 7 0C models. 





Table 9 

OGO-4 RMS of Orbital Solutions 


1 

Data Arc and Number 

of Observations 

Model 

Range 

(m) 

Range -Rate 
(cm/sec) 

Direction Cosines 
1 (mils) , m(mils) 

8/16/67 and 8/17/67 

NWL 5E-6 

84 

67 

1.7 

3.4 

R=59, RR=59, l =m=27 

SAO Ml 

63 

47 

.3 

3.5 


SAO B6.1 

60 

45 

.3 

3.4 


GSFC 1.70-C 

54 

33 

.4 

3.3 

8/17/67 and 8/18/67 

NWL 5E-6 

150 

109 

2.0 

4.0 

R=55, RR=54, ^=m=33 

SAO Ml 

90 

57 

1.4 

3.7 


SAO B6.1 

72 

49 

1.3 

3.8 


GSFC 1.70-C 

60 

28 

.9 

3.6 

8/21/67 and 8/22/67 

NWL 5E-6 

126 

89 

1.7 

3.4 

R=85, RR=82,^=m=38 

SAO Ml 

69 

51 

1.0 

3.3 


SAO B6.1 

60 

22 

.2 

1.9 


GSFC 1.70-C | 

48 

27 

.6 

2.4 

8/22/67 and 8/23/67 

NWL 5E-6 

102 

65 

.7 

3.4 

R=101, RR=98, l =m=36 

SAO Ml 

54 

51 

.7 

3.0 


SAO B6.1 

66 

39 

.7 

3.0 


GSFC 1.70-C 

i 

45 

32 

.3 

2.3 

8/26/67 and 8/27/67 

SAO Ml 

81 

57 

.8 

3.1 

R=78, RR=78, l =m=34 

GSFC 1.70-C 

72 

50 

.8 

2.9 

8/27/67 and 8/28/67 

SAO Ml 

81 

41 

.7 

3.0 

R=98, RR=97, l =m=36 

GSFC 1.70-C 

75 

50 

.8 

3.0 
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Table 10 

OGO 4 Orbital Comparisons 



Overlap Region 


8/17/67 


Model 

RMS Position Difference (meters) 

Maximum 

Radial 

Crosstrack 

Alongtrack 

Total 

Difference 

NWL 5E-6 

99.4 

407.8 

263.2 

495.4 

785.1 

SAO Ml 

31.1 

298.4 

105.7 

318.1 

484.8 

SAO B6.1 

20.4 

441.2 

131.1 

460.7 

764.9 

GSFC 1.70C 

19.9 

79.5 

84.7 

117.9 

254.7 


8/22/67 


NWL 5E-6 
SAO Ml 
SAO B6.1 
GSFC 1.70C 



151.1 

69.5 

26.7 

21.8 


151.7 
247.4 
296.2 

166.7 


215.8 

257.6 

299.0 

169.2 



8/27/67 

SAO Ml 

20.2 

25.5 

62.9 

7 

I 

.8 

144.5 


GSFC 1.70C 

4.4 

- 

8.5 

28.9 

3 

.5 

64.9 
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The best orbital fits for practically every two, day arc were obtained; when the GSFC 1.70C 
model was used. o In several 4 instances the RMS fits for the Range data we re, a, factor of two lower 
than those obtained* with, the NWL 5E.-.6 model and,; the fits to the range rate data were as much a 
factor of three lower in some cases. 

Similar results were obtained from the orbital overlap comparisons. The. results presented 
in Table* f 0; indicate thgt the smallest total, RMS position. differences were, obtained, for every case 
when, thie GSRG 1;.70-.C mpdef was used; Although, in general smaller, position differences were 
obtained, foij. all components^ thp largest reductionwas observed in the crosstrack component with 
the* GSF.O 1 .7,0£> cross, track differences on 8/22/ 67; a factor of six smaller than those obtained; using, 
the NWL 5E-6 model.. 


Although, a study of this* nature, has not been performed, for- the OGO r 6 satellite data (OGO-jB; 
orbital* paramaters are quite* similar* to those for* OGpr4) it. is anticipated thgf ; the results of such 
a study would be. comparable- to. those presented; in Tables 9. and- 10, That is^ we would expect the 
overlap errors to be reduced; significantly . 


The combined fields y?as compared to terrestrial; free -air anopialloS; and? to other geopotential, 
fields through consideration of degree variances % The degree variances computed., from the ter- 
restrial data were obtained, from, Reference 2 % The* degree yar iancee> ; cr| f,qr r the. various ; geopo- 

tential fields were, obtained, from equation (4) which is based; upon an, analysis, of; gravity appearing 
in Reference 18. Thus, 


O-J = y 2 (l - l ) 2 


2 ' - 


where 7 is the mean acceleration of gravity^ 'phpse degree* variances., fo,r* thp. input, fields^ GSRC 
1.70-C gravimetric, and* f or* the, SAO 1969 standard, earth, Reference 2p, appear in Table 1:1. 


In Table 12 there appears an RMS coefficient difference between GSFG 1.70rC and the various 
other geopotential fields. By both of these methods* of comparison GSFC 1.70rC* compares very 
well with Kohnlein T s field. G,ood agreement between GSFC 1.70-C and Kaula CA, SAO Bl#.f, 

Rapp 67, and Rapp. 68 on the, two Qieans; of comparison are obtained. 


The mean degree variance, a- 2 ; is given by equation (5) for a set of non- zonal harmonics. 
Thus, 

g l = J2 - *h, ) /2 (5) 

: T' 1 


where a^ 


log of 10 


17 



is the number of pairs of non- zonal harmonics of degree In Figure 1 a plot of the 
times the mean degree variance for* 1.70-C versus the log of the. degree is shown. 


The rule of thumb that the size of a normalized gravity coefficient of degree i is 10“ 5 /l 2 is ad- 
hered to closely by the data points in this figure. The straight line so obtained is in good agree- 
ment with one presented in Reference 19 for NWL-8 D by Anderle and Smith. 

Finally, a geoid map based upon GSFC 1.70-C and the other constants in Table 2 is given in 
Figure 2. 

Table 11 

Degree Variances o-j (m gal 2 ) 


Degree 


NWL APL 
5E6 3.5 

SAO 

Ml 

Kohnlein 

i Kaula 
656 

Kaula 

CA 

Kaula 

K8 

7.9 6.8 

7.2 

7.2 

7.5 

7.4 

7.5 

36.2 25.9 

29.0 

24.6 

26.6 

25.8 

25.5 

21.6 14.4 

16.4 

14.8 

20.6 

16.9 

17.8 

21.7 15.4 

18.2 

13.5 

19.3 

9.7 

21.1 

30.8 24.9 

17.9 

14.9 

49.6 

18.8 

19.9 

32.8 44.9 

13.3 

10.5 

6.1 


21.0 

28.4 

11.7 

8.2 

2.7 


10.4 



3.3 

3.7 




SAO Rapp 
B13.1 67 

7.4 6.8 

28.1 25.1 

19.1 14.8 

15.4 9.0 

15.1 14.9 

14.3 9.2 


Rapp 

68 


GSFC 

1.70C 


Gravimetric 


SAO 
69 S.E. 

7.4 

40.1 

19.2 
17.1 

15.7 
17.1 

6.9 

13.0 

13.8 

16.0 
11.4 

15.9 
16.1 

30.3 
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1 

GSFC 1.70-C 

and 

RMS Difference 
(unit is IQ' 6 *) 

NWL 5E-6 

.163 

APL 3„5 

.192 

SAO Ml 

.067 

Kohnlein 

.047 

Kaula .UCLA*6E>6 

.165 

Kaula CA 

.061 

Kaula K8 

.179 

SAO B6.1 

.098 

SAO B13.1 

.055 

Rapp 67 

.066 

Rapp 68 

.052 

SAO 69 S.E. 

.062 






90 



Figure 2. Geoid undulations from GSFC 1.70-C solution in meters. 


COMMENT AND CONCLUSIONS 


Several possible conclusions with regard to the performance of the combined gravity field 
for orbit computations can be made. 


1. After the proper resonant terms are inserted in the combined field, the fit to the tracking 
data is as good or better than the fits obtained with any of the fields used in composing it. 

2. The truncated field with 60% fewer coefficients fit the data almost as well as any of the 
fields tested. This implies that a 60% savings in computer time for trajectory integration could be 
realized over a (15 X 15) field with GSFC 1.70-T for some satellite prediction problems. 


3. Adjustments to the resonant coefficients were small for both GEOS- 1 and GEOS- II. This 
seems to indicate that the complement of the field with respect to these resonant coefficients 
represents the gravity field adequately. 


4. Comparisons of 1.70-C with other fields from the point of view of satellite position errors, 
degree variances, and geoid maps are favorable. 






5. The tracking data processing carried out in this study was accomplished using the NONAME 
(orbit and geodetic parameter estimation) system, Reference 21, which includes all satellite per- 
turbations such as lunar- solar .gravitation, solar radiation pressure, etc. to the required degree of 
accuracy. 

6. Professor Kaula’s conclusion concerning the superiority of the arithmetic mean of different 
solutions to any single solution has been reexamined using several additional models. 
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